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Abstract 

Particle number fluctuations and correlations in nucleus-nucleus collisions at SPS and RHIC energies 
are studied within the statistical hadron-resonance gas model in different statistical ensembles and in 
the Hadron-String-Dynamics (HSD) transport approach. Event-by-event fluctuations of the proton to 
pion and kaon to proton number ratios are calculated in the HSD model for the samples of most central 
collision events and compared with the available experimental data. The role of the experimental 
acceptance and centrality selection is discussed. 
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I. INTRODUCTION 



The study of event-by-event fluctuations in high energy nucleus-nucleus collisions opens new 
possibilities to investigate the phase transition between hadronic and partonic matter as well 
as the QCD critical point (see the reviews 1|). In our recent paper ^2^ we have presented 
a systematic study of the particle number fluctuations for pions vr = tt^ + tt^ and kaons 
K = + K~ and their correlations in the statistical model (SM) in different ensembles and 
in the Hadron- String-Dynamics (HSD) transport model Sj. In the present paper we continue 
this study and consider the fluctuations of the number of protons p= p + p, the vr-p and K-p 
correlations, and fluctuations of the p/vr and K/p particle number ratios. The HSD results 
for the fluctuations in the particle number ratios are compared with NA49 data for Pb+Pb 
collisions at SPS energies ^] for the proton to pion ratio as well as preliminary STAR data for 
Au+Au collisions at RHIC energies for the proton to pion and kaon to proton [6] ratios. 

The paper is organized as follows. In Section II the observables for particle number fluctu- 
ations and correlations are introduced. In Section III the SM and HSD results for the scaled 
variances and correlation coefficients of particle number fluctuations are presented for cen- 
tral nucleus-nucleus collisions at SPS and RHIC energies. In Section IV the HSD transport 
model results for the fluctuations of the p/n and K/p ratios are presented and compared to 
the available data on proton to pion ratio fluctuations. The role of the centrality selection and 
experimental acceptance is discussed using the HSD results. A summary closes the paper in 
Section V. 

II. OBSERVABLES FOR PARTICLE RATIO FLUCTUATIONS 

We define the covariance for particle species A and B as: 

A{Na,Nb) = (AN^ANb) = (NaNb) - {Na){Nb) , (1) 

where (Na) is the event-by-event average of the particle multiplicity Na and AA^^ = Na — {Na)- 
The scaled variance uja and correlation coefficient pab are defined by: 

_ A{Na,Na) _ {ANaANb) 

= TT^TT , PAB = ; — 7^ ■ l^J 
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2 _ {{ARabY) 



(4) 



Furthermore, the fluctuations in the ratio Rab = Na/Nb will be characterised by: 

JrIb)' 

One finds to second order in ANa/ {Na) and ANb/ (Nb): 

r 1 1/2 

2 ^ ^A ^B r) ^a^b 

" {Na) + {Nb) ^^^^ [{Na){Nb)_ 

We recall that the experimental data for Na/Nb fiuctuations are usually presented in terms 
of the so called dynamical fiuctuations l?!^ 

dyn ■ f 2 2 \ \ 2 2 1 1/2 /r\ 

^AB = ^^S^y'^AB ~ ^AB,m,ix)\^AB ~~ ^AB,mix\ ' \^) 

where cr^^ is defined by Eq. (jlj) while cr^^ corresponds to the mixed events background (see 
Ref. fl): 

2 _ 1 1 

^AB,rmx " ^ + ■ (6) 

III. PARTICLE NUMBER FLUCTUATIONS AND CORRELATIONS 

In this section we present the results of the SM and HSD for the particle number fiuctuations 
and correlations in central nucleus-nucleus collisions. The procedure of the calculations is 
essentially the same as in our previous paper j^. The SM and HSD results for scaled variances 
of kaons K = + K~ , pions tt = 7r+ + vr", and protons p= p + p, and the correlation 
coefficients pKn, Pnp, and pxp (defined by Eq. ([2])) are presented in Tables I and II. Figures [H 
and [2] show ujp and Pp^, pxp, respectively. We note the results for ujk, ^-k and pk-k have been 
presented in Ref. Q]. 

Let us first comment on the SM results for the particle number fiuctuations and correlations 
in different ensembles. The global charge conservation laws suppress the particle number fiuc- 
tuations in the canonical ensemble (CE) in a comparison with the grand canonical ensemble 
(GCE). Exact energy conservation in the micro canonical ensemble (MCE) makes this sup- 
pression even stronger. This leads to the result that the MCE scaled variances uja become 
smaller than unity. The contributions to the correlation parameter pab in the SM stem from 



1 Other dynamical measures, such as $ @, 0| and F [13], may be also used. 
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TABLE I: The chemical freeze-out parameters T and for central Pb-|-Pb (Au-|-Au) collisions along 



the chemical freeze-out line The hadron-resonance gas model results are presented for the number 
density olp+p, (particle number densities are the same in all statistical ensembles in the large volume 
limit), scaled variance Wp, and correlation parameter pp.,^ and pKp in the GCE, CE, and MCE. The 
SM results for 71,^, nx, Wtt, and px-w are presented in Table I of Ref. 



two sources: resonance decays and global conservation laws. Resonances decaying into pairs of 
particles of species A and B produce a corresponding correlation and yield a positive contribu- 
tion to pab- On the other hand, the conservation laws lead to an anticorrelation and thus to 
negative contributions to pab- 

As seen from Figs. [Hand [2] the HSD results for ojp, pp^, and pxp (denoted by the solid lines) 
are close to the CE and MCE results for low SPS energies. One may conclude that the influence 
of conservation laws is more stringent at low collision energies. The same conclusion follows 
from the HSD results for uj^^, ujk and pk-k (see Figs. 1 and 2 in Ref. 2[]). The HSD values for 
uja and pab increase, however, at high collision energies and a sizeable deviation of the HSD 
results from those in the MCE SM is observed with increasing energies for y^J/vw > 200 GeV. 

We point out again that important aspects of the event-by-event fluctuations in nucleus- 
nucleus collisions are the dependence on the centrality selection and experimental acceptance. 
We accordingly discuss the role of these effects using HSD results for the scaled variance uj.„ of 
the pion number fluctuations. In Fig. [3] the scaled variance cUtt (calculated within HSD) for the 



4 





HSD full acceptance 


[GeV] 




{Nk) 


{Np) 




UK 


UJp 


PKit 


Ppn 


PKp 


6.27 


612.03 


43.329 


181.83 


0.961 


1.107 


0.506 


-0.091 


0.048 


-0.137 


7.62 


732.11 


60.801 


180.03 


1.077 


1.141 


0.526 


-0.063 


0.025 


-0.128 


8.77 


823.71 


75.133 


179.43 


1.159 


1.168 


0.546 


-0.033 


0.016 


-0.129 


12.3 


1072.3 


116.44 


180.81 


1.378 


1.250 


0.596 


0.046 


-0.006 


-0.126 


17.3 


1364.6 


165.52 


186.97 


1.619 


1.348 


0.641 


0.126 


-0.010 


-0.118 


62.4 


2933.9 


449.29 


240.53 


3.006 


1.891 


0.863 


0.412 


0.074 


-0.029 


130 


4304.2 


692.59 


307.31 


4.538 


2.378 


1.020 


0.557 


0.177 


0.067 


200 


5204.0 


861.77 


352.91 


5.838 


2.765 


1.122 


0.634 


0.251 


0.135 



TABLE II: The HSD results for the average multiphcities (Nt^), {Nk), {Np) and values of uj,^, lok, 
and pKn for central (impact parameter 6 = 0) Pb-|-Pb (Au-|-Au) collisions at different cm. energies 

y/SNN- 



full acceptance and for the experimental acceptance are shown in nucleus-nucleus collisions for 
zero impact parameter 6 = (see the next Section for details of the experimental acceptance). 
Introducing the probability q of pion experimental acceptance as the ratio of an average accepted 
to the total multiphcities, q = {N^^'^) / {N^"^) , one finds: 

a;^ = 1 - q + q cuf " . (7) 

Equation ([7]) connects the scaled variance cj^"'' in the full An space with cj"^'^ defined for the 



experimental acceptance. The acceptance scaling ([7]) assumes (see, e.g. Ref. |l2j) an absence 
of particle correlations in momentum space. Figure [3] demonstrates that the acceptance scaling 
([7]) underestimates the scaled variance u^'^'^ at RHIC energies. 

The samples of collision events selected experimentally are 3.5% of most central collision 
events in Pb-|-Pb collisions at the SPS energies and 5% in Au-|-Au collisions at RHIC energies. 
Figure m presents the HSD results for u;^ in these samples of the most central events and their 
comparison with the HSD results at zero impact parameter. One finds much larger values of 
ti^Tr in the centrality selected samples than for 6 = 0. The effect is especially strong at RHIC 
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FIG. 1: (Color online) The SM results in the GCE, CE, and MCE ensembles and the HSD results 
(impact parameter 6 = 0) are presented for the scaled variance ujp in Pb+Pb (Au+Au) collisions at 
different cm. energies ^snn- 

energies. This can be qualitatively understood within the model of independent sources: 

= + Up , (8) 

where is the scaled variance for pions produced by one source, ujp is the scaled variance for 
the fluctuations of nucleon participants, and is the average pion multiplicity per participat- 
ing nucleon which increases monotonously with collision energy. Collisions with zero impact 
parameter correspond to up = 0. Thus, can be approximately taken as aX b = 0. The 
HSD results correspond approximately to Up = 0.5 for the 3.5% most central Pb+Pb collisions 
at SPS energies and up = 1 for the 5% most central Au+Au collisions at RHIC energies. Please 
note that we used the restrictions on impact parameter b in the HSD calculations to form the 
samples of most central events. The results of the model of independent sources ([8]) for are 
shown by the dashed line in Fig. H] and are close to the actual values of the HSD simulations 
for cUtt. 
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FIG. 2: (Color online) The SM results in the GCE, CE, and MCE ensembles and the HSD results 
(impact parameter 6 = 0) are presented for the correlation parameters ppT^ {upper pannel) and pxp 
{lower panel) in Pb+Pb (Au+Au) collisions at different cm. energies ^/sNN■ 
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FIG. 3: (Color online) The HSD results for lo-,^ are presented for Pb+Pb (Au+Au) collisions with 
zero impact parameter (6 = 0) at different cm. energies ^snn- The upper solid line corresponds 
to the full 47r-acceptance and the middle one to the experimental acceptance. The lower dashed line 
corresponds to the acceptance scaling Eq. ([7]). 

IV. FLUCTUATIONS OF PARTICLE NUMBER RATIOS 

A comparison of the SM results for fluctuations in different ensembles with the data looks 
problematic at present; the same is true for most other theoretical models. This is because of 
difficulties in implementing the experimental acceptance and centrality selection which, how- 
ever, can be taken into account in the transport approach. In order to compare the HSD 
calculations with the measured data, the experimental cuts are applied for the simulated set 
of HSD events. In Fig. [5] the HSD results of a"^^^ ([5]) for the p/rc and K/p ratios are shown in 
comparison with the experimental data by the NA49 Collaboration at the SPS [4] and the pre- 
liminary data of the STAR Collaboration at RHIC The results of the UrQMD calculations 
for a^Y^ at the SPS energies (from Ref. \^\) are also shown by the dashed line. 



For the SPS energies we use the NA49 acceptance tables from Ref. j4|. For the RHIC 
energies we use the following cuts: in pseudorapidity, \ri\ < 1, and in transverse momentum. 
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FIG. 4: (Color online) The HSD results for uj.,^ for Pb+Pb (Au+Au) collisions at different cm. energies 



and the upper one to the experimentally selected samples of collision events. The dashed line reflects 
the model of independent sources ([8]). 



0.2 < pt < 0.6 GeV/c for kaons and pions and 0.4 < px < I GeV/c for protons [5|, [g, Il3| . 
We note also, that HSD results presented in Fig. [S] correspond to the centrality selection as 
in the experiment: the NA49 data correspond to the 3.5% most central collisions selected via 
veto calorimeter, whereas in the STAR experiment the 5% most central events with the highest 
multiplicities in the pseudorapidity range |?7| < 0.5 have been selected. 

At the SPS energies the HSD simulations lead to negative values of a'^^^ for the proton 
to pion ratio. This is in agreement with the NA49 data in Pb+Pb collisions. On the other 
hand HSD gives large positive values of 0"^!^" at RHIC energies which strongly overestimate 
the preliminary STAR data for Au+Au collisions [5]. For a'j^^ only preliminary STAR data in 
Au+Au collisions are available which demonstrate a qualitative agreement with the HSD 
results (Fig. [5]). The HSD results for a'j^p show a weak energy dependence in both SPS and 
RHIC energy regions. A peculiar feature is, however, a strong 'jump' between the SPS and 
RHIC values seen in the lower panel of Fig. [5], in the HSD calculations which is caused by the 



^SNN within the full 47r-acceptance. The lower solid line corresponds to zero impact parameter (b = 0) 
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HSD: (p + p) / {n + n) 
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FIG. 5: (Color online) The solid lines present the HSD results of cr'^^" (jS]) for the (p + p)/{'^~^ + vi" ) 
{upper panel) and {K + K^)/{p + p) {lower panel) ratios. The data are from Refs. [J,l5|,|6|. 



different acceptances in the SPS and RHIC measurements. 

The influence of the experimental acceptance is clearly seen at 160 A GeV where a switch 
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from the NA49 to the STAR acceptance leads to the jump in cx^^ by 3% - lower panel of Fig. [5l 
On the other hand, our calculations for Pb+Pb (3.5% central) and for Au+Au (5% central) 
collisions - performed within the NA49 acceptance for both cases at 160 A GeV - shows a very 
week sensitivity of a'f^^ on the actual choice of the collision system and centrality - cf. the 
coincident open circle and triangle at 160 A GeV in the lower panel of Fig. [5l 



V. SUMMARY AND CONCLUSIONS 



We have studied the event-by-event fluctuations of the number of protons (and anti-protons), 
the proton-pion and proton-kaon correlations, and the fluctuations of proton to pion and kaon 
to proton ratios in central Pb-|-Pb and Au-|-Au collisions from low SPS up to top RHIC energies. 
The analysis has been performed within the statistical hadron-resonance gas model for different 
statistical ensembles - the grand canonical ensemble (GCE), canonical ensemble (CE), and 
micro-canonical ensemble (MCE) - and in the Hadron-String-Dynamics transport approach. 
We have found that the HSD results at SPS energies are close to those in the CE and MCE 
statistical model. This indicates a dominant role of resonance decays and global conservation 
laws at low energy nucleus-nucleus collisions. On the other hand, substantial differences in HSD 
and statistical model results have been observed at RHIC energies which can be attributed to 
non-equilibrium dynamical effects in the HSD simulations. 

The HSD results for a^^" appear to be close to the NA49 data at the SPS. The data for a'^p 
in Pb+Pb collisions at the SPS energies will be available soon and allow for further insight. 
A comparison of the HSD results with preliminary STAR data in Au+Au collisions at RHIC 
energies are not fully conclusive: cr'^^" from HSD calculations are approximately in agreement 
with data for kaon to proton ratio, but overestimate the experimental results for proton 
to pion ratio. New data on event-by-event fluctuations in Au+Au at RHIC energies will help 
to clarify the situation. 
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